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Abstract 

Recent advances in N-body simulations of cold dark matter halos 
point to a substantial density enhancement near the center. This means 
that, e.g., the 7 ray signals from neutralino dark matter annihilations 
would be significantly enhanced compared to old estimates based on 
an isothermal sphere model with large core radius. 

Another important development concerns new detectors, both space- 
and ground-based, which will cover the window between 50 and 300 
GeV where presently no cosmic 7-ray data are available. 

Thirdly, new calculations of the 7-ray line signal (a sharp spike of 
10~^ relative width) from neutralino annihilations have revealed a hith- 
erto neglected contribution which, for heavy higgsino-like neutralinos, 
gives an annihilation rate an order of magnitude larger than previously 
predicted. 

We make a detailed phenomenological study of the possible detec- 
tion rates given these three pieces of new information. We show that 
the proposed upgrade of the Whipple telescope will make it sensitive 
to a region of parameter space, with substantial improvements possible 
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with the planned new generation of Air Cherenkov Telescope Arrays. 

We also comment on the potential of the GLAST satellite detector. 

An evaluation of the continuum 7-rays produced in neutralino an- 
nihilations into the main modes is also done. It is shown that a com- 
bination of high-rate models and very peaked halo models are already 
severely constrained by existing data. 



1 Introduction 



The dark matter problem is one of the most outstanding problems con- 
fronting cosmology and astrophysics today. The question of the nature of 
the dark matter in the Universe is a cross-disciplinary one which may need 
elements of particle physics for its solution. With new observations, the 
possible candidates get more and more constrained. 

Let us first recall that from the particle physics point of view, the theo- 
retically preferred (Einstein-De Sitter) Universe has the simple description 

where as usual we have introduced the normalization of the energy density 
to the critical density 

Peru 1.9- 10-29/^2 gcni-3' 

and SIa is the contribution to the energy density from a cosmological con- 
stant (or equivalently from vacuum energy). 

The Einstein-De Sitter model (||) has the attractive features that it is 
simple, avoids finetuning, and may be explained by a period of inflation in 
the earliest Universe. Since Big Bang nucleosynthesis (BBN) puts an upper 
limit to the baryonic contribution of 

Qfc/i^ < 0.026 , (3) 

with h related to the Hubble constant Hq hy h = Hq/{100 kms~^Mpc^^) 
(observationally, h lies between 0.4 and 0.8), non-baryonic dark matter dom- 
inates the energy density by a large factor in this type of model. Indeed, 
there is a number of different observations on scales from dwarf galaxies and 
upwards which point to a larger value of Qtot than that allowed by nucle- 
osynthesis [^. It therefore seems that non-baryonic dark matter has to be 
present in substantial amounts in the Universe today. 

One of the prime candidates for the non-baryonic component is provided 
by the lightest supersymmetric particle, plausibly the lightest neutralino 
X Q. Supersymmetry seems to be a necessity in superstring theory (or M- 
theory) which unites all the fundamental forces of nature, including gravity. 
In most versions of the low-energy theory there is a conserved multiplicative 
quantum number, R-parity, which makes the lightest supersymmetric parti- 
cle stable. Thus, pair-produced neutralinos in the early Universe which left 
thermal equilibrium as the Universe expanded should have a non-zero relic 
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abundance today. If the scale of super symmetry breaking is related to that 
of electroweak breaking, may come out in the right order of magnitude 
to explain the non-baryonic dark matter. 

In addition, neutralinos are generically found to decouple at a temper- 
ature that is roughly M^/20, which means that they are non-relativistic 
already at decoupling and certainly behave as Cold Dark Matter (CDM) by 
the time of matter dominance and structure formation. Most analyses of 
large scale structure formation in the Universe indeed find that CDM mod- 
els (perhaps with a small addition of, e.g., massive neutrinos) give the best 
description of observational data. 

Although accelerators like the Tevatron at Fermilab and LEP at CERN 
have already started to probe regions of parameter space, so far without 
finding any super symmetric signals, the constraints imposed are not yet 
very restrictive. Adding the requirement that neutralinos make up most of 
the dark matter, one finds that there are viable models in the mass inter- 
val 30 GeV < < 10 TeV. In the low mass range, the direct detection 
at accelerators or in terrestrial detectors sensitive to the weak interactions 
generated by x particles of the Milky Way halo as they pass the Earth are 
probably the most promising methods of detection (see |Q for an extensive 
review) . 

At the high mass end (or, depending on halo paramaters, maybe even 
over the full mass range) indirect detection methods are competitive. The 
best indirect signals are given by neutrinos from the Sun or the central 
region of the Earth, and by almost monoenergetic 7 rays lines originating 
from neutralino annihilations in the Milky Way halo. 

Neutrinos can escape from the centre of the Sun or Earth, where neu- 
tralinos may have been gravitationally trapped and therefore their den- 
sity enhanced. Gamma rays may result from loop-induced annihilations 
XX 77 [i or XX ^ |]. 

The rates of these processes are difficult to estimate because of uncer- 
tainties in the supersymmetric parameters, cross sections and halo density 
profile. However, in contrast to other proposed detection methods they have 
the virtue of giving very distinct, "smoking gun" signals: high-energy neu- 
trinos from the centre of the Earth or Sun, or monoenergetic photons with 

= My. or = M^(l — r?T,^/4M^) from the halo. The neutrino signal has 
been thoroughly discussed in the literature [^. In this article, we concen- 
trate on 7 ray lines, where three recent lines of development have prompted 
us to take a new look at this process (for early studies, see 0). 

Firstly, a number of numerical studies using N-body simulations ^, 
10, 11 1 have shown that in CDM cosmologies the hierarchical way that halos 
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form means that they should have a density profile which is quite steep 
with a singular behaviour ~ r~" with estimates of a ranging from 0.2 to 
1.7. The physics underlying this behaviour is simple: in CDM models with 
scale-invariant gaussian initial fluctuations small clumps go nonlinear first, 
where the large average background density at the early epochs means that 
the smallest clumps are densest. Through the amplifying action of gravity, 
larger regions of overdensity form by the merging of smaller substructures. 
Through tidal interactions the overall halo distribution gets smeared out, 
but near the center some of the densest clumps leave a trace in form of 
their high average overdensity. Taking the interplay between the CDM and 
baryonic parts of the halo into account, a satisfactory description of the 
rotation curves of even dwarf spheroidal galaxies can be achieved 

Since the probability for two CDM particles such as neutralinos to meet 
and annihilate each other is proportional to the square of the neutralino 
density, this enhancement near the galactic center could be of utmost im- 
portance for the 7 ray line signal. 

The second development concerns the calculation of the processes XX ~^ 77 
and XX ~^ -^7) which only recently were computed fully to one loop in the 
MSSM, the minimal supersymmetric extension of the Standard Model of 
particle physics [@, In these calculations it was shown that previously 
neglected contributions could increase the predicted rates by an order of 
magnitude, especially at the high-mass end of the allowed range. 

In the high mass range, > 250 GeV, existing Air Cherenkov Tele- 
scopes (ACT) could already be sensitive to the 7 ray lines. However, the 
third interesting development is the planning and construction of a new 
generation of ACTs which will have larger area, lower threshold and better 
energy resolution. This will improve substantially the discovery potential 
of these telescopes for supersymmetric dark matter. Also, plans are being 
made for a satellite-borne detector (GLAST is one of the design concepts) 
which can compensate a factor ~ 10^ smaller area than ACTs by a larger 
angular acceptance, better energy resolution, lower energy threshold and 
longer integration time. 

The plan of this article is as follows. In Section ^, we give an introduction 
to the MSSM and the predictions of the 7 line flux. In Section ^ we discuss 
various halo models and their implications for the 7 ray flux in ACTs from 
neutralino annihilations. In Section |^ we give an overview of present and 
planned detectors for cosmic 7 rays, and in Section || we discuss the detection 
potential of several existing or upcoming ACTs. In Section ^ we discuss 
the discovery potential of satellite detectors, and in Section |^ we comment 
also on the diffuse continuum 7 rays that would be produced as secondary 
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annihilation products if neutralinos make up a significant part of the dark 
matter. FinaUy, in Section |^ we discuss our results and give some concluding 
remarks. 

2 Overview of MSSM Results for 7 Lines 

We have performed a detailed phenomenological analysis of the annihilation 
rate of non relativistic neutralinos into two photons and into a photon and 
a Z boson. The two processes, which have the experimentally significant 
feature of giving nearly monochromatic photons in the final state, with an 
energy ~ and ~ ~ -^I/^M^), respectively, have been 

studied in Refs. |^] and |5| where for the first time a full one loop calculation 
of the two cross sections was performed. Considering a broad selection of 
MSSM models in which the neutralino is the lightest supersymmetric particle 
and a good dark matter candidate, the aim of this Section is to examine in 
which cases the two annihilation processes are important. We will compare 
the relative values of the cross sections in view of extracting information on 
the nature of the neutralino from a possible detection of one or both the 
gamma ray lines from neutralino annihilations in the galactic halo. 

In the minimal = 1 supersymmetric extension of the Standard Model 
four neutral spin-1/2 Majorana particles are introduced, the partners of 
the neutral gauge bosons i?, and the neutral CP-even higgsinos Wl, 
Diagonalizing the corresponding mass matrix, four mass eigenstates 
are obtained: 

X° = NaB + Ni2W^ + NisH^ + Ni^H° (4) 

The lightest of these, Xi or simply x, is commonly referred as the neutralino. 
It is useful to introduce the gaugino fraction Zg defined as 

Zg = \Nn\^ + \Ni2\^ (5) 

and classify the neutralino as higgsino-like when Zg < 0.01, mixed when 
0.01 <Zg< 0.99 and gaugino like if Zg > 0.99. 

The general R-parity conserving version of the MSSM is defined by 63 
free parameters |^]. Some simplifying assumptions are necessary to reduce 
the number of parameters and get to a format which is numerically tractable 
(for a complete discussion of this procedure see Ref. m, 1). In the 
scheme we use there are seven free parameters: the higgsino mass parameter 
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1000 GeV < 1^1 < 30000 GeV 
1000 GeV < I Ma I < 50000 GeV 
1.2 < tan/3 < 50 
OGeV < niA < 10000 GeV 
100 GeV < mo < 30000 GeV 
—3 mo < At, < 3 mo 
—3 mo < At < 3 mo 



Table 1: Considered range for the free parameters in our realization of the 
Minimal Supersymmetric Standard Model. 

fi, the gaugino mass parameteiQ M2, the ratio tan (3 of the vacuum expecta- 
tion values of the two Higgs fields H2 and , the mass of the pseudoscalar 
Higgs ruA, the sfermion mass parameter mo and other two soft supersym- 
metry breaking parameters Af, and Af. Our analysis is based on 16 scans 
over this parameter space. Some of them are broad samplings in the param- 
eters, others are dedicated to more specific portions of the parameter space 
that are interesting for the present purpose. Table |l] contains the maximal 
and minimal values for each parameter; in the given intervals, or in most 
cases in subintervals, /i and M2 are scanned logarithmically, tan/3 is mostly 
scanned linearly and in a few scans logarithmically, while all the other pa- 
rameters are varied along a linear scale. Rejecting those models that violate 
all known experimental bounds, we are left with a sample of about 62,500 
supersymmetric models which could describe Nature. For each of these, the 
neutralino relic density Q^h"^ has been computed using the exact treatment 



of resonances |15] and including coannihilations between neutralinos and 



charginos [16|. We restrict to those models which can give the neutralino 



as the main constituent of galactic halos, selecting the models which have 
a relic density that is in the interval 0.025 < fi^/i^ < 1 (the upper limit is 
given by the condition of not overdosing the universe, or equivalently giving 
a too low age of the Universe). About 22,000 models in our sample fulfill 
this condition and this is the selected sample that is shown in the graphs in 
this paper. 

In Fig. m and Fig. |2| we show the predictions for the annihilation rate of 
neutralinos into two photons and into a photon and a Z boson as a function 
of the neutralino mass, while in Fig. |3|they are shown one versus the other. 

^The usual Grand Unification Theory (GUT) relations for gaugino masses are assumed 
in all the results presented, their release is not expected to give major modifications 
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It is useful to compare 2 va2'y and vaz-y as the 7 line flux is proportional to 
these quantities (the 2 comes from the fact that there are two photons in that 
final state) . In these figures as in all those that follow, each marker represents 
a model in our sample, we have however arbitrarily fixed a maximal density 
of points which are printed (it might be worth reminding that the point 
density in these kind of figures is meaningless and just depends on the way 
the parameter space is sampled). The visual effect is the same as drawing 
shaded zones but in this way we avoid the problem of deciding whether or 
not to include isolated points or empty sub-regions. 

A common feature of 2va2'y and vaz'y is that the highest values are 
given by very pure higgsino-like neutralinos (we have included models with 
Zg as small as 10~^), with a mass around 500 GeV. As it was already shown 
in Refs. @ and |^, the cross sections tend to a constant value of about 
10"^*^ cm^s~^ in the limit Zg ^ and » Mw- The higgsino-like 
neutralino branch starts at very high masses, where the relic density is close 
to 1 (coannihilations of neutralinos with the lightest chargino are crucial to 
show that neutralinos with masses highier than 3 TeV can be cosmologically 
acceptable |jl^ ) , and continues with increasing cross sections down to about 
500 GeV where f^^/i^ approaches 0.025. This is mainly due to the large 
cross section for the annihilation into a W boson pair. In fact, when this 
process becomes kinematically forbidden (for neutralino masses below the 
mass of the W boson) higgsino-like neutralinos become again cosmologically 
interesting candidates and, together with mixed neutralinos, again give the 
highest cross sections in that mass region. 

Whereas for higgsino-like neutralinos 2 va2'y is in most cases greater than 
vaz^y, it is generally the opposite for mixed neutralinos. For this class of 
neutralinos in the mass interval 100 — 500 GeV the branching ratios into 
the Z7 final state is in most cases dominant over twice the branching ratio 
into 27. Another interesting feature is that over the whole mass scale the 
value of vaz'y for mixed neutralinos is concentrated in a band not wider than 
two orders of magnitude, while the values of va2'y are more homogeneously 
distributed over 5 to 6 decades. This is visible also in Fig. ^ where the mixed 
neutralino band starts at about 6 — 8-10"^^ cm^s~^ for both 2 va2y and vaz^y 
and extends down to to a value of 2va2'y of lO^^^cm^s^^ for models that 
can give a vaz^y as high as 10~^^ cm^s~^. We illustrate this point in Fig. ^ 
where we consider a scan of models with mixed neutralinos of masses in the 
range 150 — 400 GeV. In these figures we have indicated which is the class 
of diagrams that gives the main contribution to the cross sections (all the 
diagrams involved in the computation are shown in Refs. and Q). In 
the case of vazy almost always the W boson-chargino loop diagrams are 
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dominant over the other diagrams. For 2 f C727 instead, the contributions for 
the fermion-sfermion loop diagrams can be more frequently of the same order 
of magnitude as the W boson-chargino diagrams and severe cancellations 
can take place, with usually a small cross section when the fermion-sfermion 
diagrams are dominant. From Fig. ^ it is also evident that Higgs-chargino 
loop diagrams and Goldstone boson-chargino loop diagrams do not play a 
main role and this is true for all classes of neutralinos over the whole mass 
range. 

Given the large size of our sample of scanned models, we expect the 
narrowness of the band for large neutralino masses to be a real effect 
which could increase the discovery potential of present and future ACTs. 
In the heavy mass range, it may be appropriate to sum the 77 and 
contributions, since to resolve the two lines an energy resolution better than 
around 1 % would be needed. 

Gaugino-like neutralinos give cross sections that rarely exceed, for both 
processes, a few times 10"^*^ cm^s~^. In this case the cross sections are 
poorly correlated and the large spread in Fig. |^ is due mainly to gaugino- 
like neutralinos. 

A final comment can be made about low mass neutralinos. Here vaz^-j 
goes to zero at the threshold M-^ Mz/2 as expected because this limit 
corresponds to having no photon in the final state and the annihilation of 
neutralino pairs into a Z boson is forbidden by C parity (for clarity we have 
excluded models with neutralino masses lower than 70 GeV from Fig. ^). 
In addition, 2 va2'y has low values as well because the contribution from W 
boson-chargino loop diagrams becomes small (there is no contribution from 
the imaginary part and the real part is generally large only for heavy pure 
higgsino-like neutralinos) . 

3 Milky Way halo models 

The photon flux from neutralino annihilation in the galactic halo critically 
depends on the neutralino distribution. Since the flux per unit volume is 
proportional to the neutralino density squared, any enhancement of the den- 
sity would result in a sharp increase in the photon flux; such an enhancement 
is provided by dark matter radial halo profiles that are peaked towards the 
galactic center. In this section we will examine in detail this possibility and 
compute the full dependence of the flux on the galactic coordinates in some 
selected cases. 

The mass distribution in the Milky Way and the relative importance 
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of its three components, the bulge, the disk and the halo, are poorly con- 
strained by available observational data. Although the dynamics of the 
satellites of the galaxy clearly indicates the presence of a non-luminous mat- 
ter component, a discrimination among the different radial dark matter halo 
profiles proposed in the literature is not possible at the time being |17|. Our 



approach is to assume that dark matter profiles are of a universal functional 
form and to infer the Milky Way dark matter distribution from the results of 
N-body simulations of hierarchical clustering in cold dark matter cosmolo- 
gies. The predicted profiles in these scenarios have been tested to a sample 
of dark matter dominated dwarf and low-surface brightness galaxies which 
provide the best opportunities to test the spatial distribution of dark matter. 
Actually this field of research is in rapid evolution and slightly discrepant 
results have recently been presented fs], 11]; we will concentrate on the 



Kravtsov et al. profile [10| and the Navarro et al profile Q, comparing the 
results with the canonical modified isothermal sphere profile which is widely 
used in dark matter calculations. 

Among the general family of dark matter halo profiles 

^^"■^ (r/a)7 [l + (r/a)«](/3-7)/-' 

Kravtsov et al. propose a coreless profile with a mild singularity towards 
the galactic center with 7 ~ 0.2—0.4; taking into account these two extremes, 
we will refer to the Ka profile as to the one which is defined by (a, f3, 7) = 
(2, 3, 0.2), while the Kb profile is chosen as (a, (3, 7) = (2, 3, 0.4). In previous 
work 1^], Navarro, Frenk and White had obtained a cuspy profile which 
has (a,/3, 7) = (1,3,1) (hereafter NFW profile). Both results are not in 
agreement with the modified isothermal distribution, (a, /3, 7) = (2, 2, 0) 
(hereafter Sp profile), extensively used in the analyses of observed rotation 
curves, which is non singular with a core of radius ~ a. 

There are models [^8| which have a more singular behavior near the 
galactic center, and which would give enormously enhanced rate in that 
direction. However, there is observational evidence against such singularities 
from cluster gravitational lensing and the rotation curves of dwarf spiral 
galaxies ||T^. On the other hand, the discrepancy between the 1/r central 
cusp of the NFW profile and the experimental data from the dwarf spheroidal 
DDO 154 has been explained in Ref. assuming an additional component 
of dark baryons. Of course, also a moderately steeper profile like 1/r^'^, 
which is found in some simulations and which has some theoretical 
underpinning [^], may be acceptable. That would give higher estimates 
than the one we use for the neutralino-induced fluxes, but we choose to be 
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more conservative. 

We will consider only the case of spherical profiles; introducing a flat- 
tening parameter may enhance the value of the flux but the effect is not 
expected to be dramatic for this neutralino detection method and we prefer 
not to introduce another factor of uncertainty. 

The normalization constant of the halo profile, which we choose as the 
value of the halo density po at our galactocentric distance Rq, the core 
radius a and Rq itself are the three parameters which will be relevant to 
give a prediction for the gamma line flux. N-body simulations indicate that 
there may be a correlation between normalization constant and core radius, 
but we are not in the position of extracting a firm numerical prediction 
for the Milky Way. After choosing the functional form of the halo profile, 
we determine the regions in the parameter space which are allowed by the 
assumptions that follow below. 



We follow Dehnen and Binney |17] and assume that the total mass of 
the galaxy inside 100 kpc is 7 it 2.5 x lO^^M©. This estimate is obtained by 
combining the result of Kochanek [^] who used the observational constraints 
from the velocity distribution of the Milky Way's satellites, the local escape 
velocity of stars and the Local Group timing model, and the result inferred 
in Ref. [^] from the dynamics of the Magellanic Clouds and Stream. Both 
estimates are to some extent model dependent and this is the reason for the 
large error bound chosen. We then suppose that, for each matter distribution 
model, on average ten per cent of the total mass inside 100 kpc is due to 
the contributions of the disk and the bulge (value obtained from Table 4 in 
[p!?! ) and derive a value for the the mass of dark matter halo inside 100 kpc 

Mh{r < 100 kpc) = (6.3 ± 2.5) x IO^^Mq (7) 

A second constraint is given in terms of the contribution of the halo to 
the local rotation curves which can be determined in terms of the observed 
value of the local rotation velocity Qq and from the measurements of the 
contribution of the disk (which is dominant over the contribution from the 
bulge) : 

vliRo) = el-vj{Ro) (8) 

The lAU standard values for Qq and our galactocentric distance Rq, based 
on a review by Kerr and Lynden-Bell in 1986, are respectively 220 it 
20kms~^ and 8.5 it 1.1 kpc. The two quantities are not independent vari- 
ables, their ratio is fixed by local stellar kinematic measurements of the 
Oort constant. In more recent estimates lower values have been obtained: 
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Reid m has found Rq = 7.2 ± 0.7 kpc, while in a very recent work Oiling 
and Merrifield [||] has found Rq = 7.1 ± 0.4 kpc and Go = 184 ± 8kms-^ 
Even the value of the contribution of the disk to the local rotation velocity 
has large uncertainties and is to some extent model dependent. The value 
we can infer from Kuijken and Gilmore |27|, and Gould [28|, who examine 



a data set of velocities of local stars as a function of their height over the 
galactic plane, is in the range Vd{Ro) = 137 it 16kms~^. Maximum disk 
models can give a higher Vd{Ro), while in thin disk models the value can 
be lower. We will accept the quoted value without considering disk models 
in detail. Combining the values of @q and Vfi{Ro) we obtain the following 
ranges of values of Vh{RQ) for a given galactocentric distance Rq: 

VhiRo) ~ 128 - 207kms~^ if Rq = 8.5 kpc (9) 

VhiRo) ~ 86- 149kms"^ if iio = 7.1kpc (10) 

In Fig. I we show, for each of the halo models considered, the regions in 
the plane (a, po) which are compatible with the limits given in Eq. ^ and 
Eq. (P) , ( |To|) . A common feature of all the profiles is that the allowed regions 
are shifted to lower values of the core radius increasing the galactocentric 
distance from 7.1 kpc to 8.5 kpc. Focusing for instance on the Ka profile 
with Rq = 8.5 kpc, the values of a and po which are compatible with our 
constraints fall in an area with a diamond-like shape, where the upper left 
bound and the lower right bound are given respectively by the maximum and 
minimum values in Eq. (|^ , while the upper right and lower left are due to 
Eq. d^. A similar behavior is seen in all the other cases. For the Sp profile a 
core radius equal to zero is not excluded, but this is an extremely unrealistic 
case and we have decided to arbitrarily fix a lower value of a = 1 kpc. 

The gamma ray flux from neutralino annihilation in the galactic halo is 
given by: 

$7W = fl^/ p\l)dlW (11) 

47rM^ J line of sight 

where ip is the angle between the direction of the galactic center and that 
of observation, and where = 2 for ^ 77j -^7 = 1 for xx ~^ -^^T- 
Separating the dependence on the halo model from the part which is related 
to the the values of the cross section and the neutralino mass, we rewrite it 
as 

,-11 f ^7 



^^(tp) ~ 1.87-10"" ^ — ^— 

V 10-29 cm^s-i, 

- J(V') cm"2 s"^ sr"^ (12) 
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where we have defined the dimensionless function 

J (4,) = — — — / p\i) dim . (13) 

8.5kpc \0.3GeY /envy June of sight ^ ' ^ ' 

Picking in Fig. ^, for the four halo models, a couple of values (a, po) which 
gives rather conservative results and choosing Rq = 8.5 kpc, we show in 
Fig. |6| in each case the angular dependence of the function J. The profiles 
have been considered valid up to the capture radius of the black hole at the 
galactic center which is about 0.01 pc for a mass of Sgr A* M ~ 10® Mq. 



we plot ( J (0) ) versus 
and for both choices of 



Modifications to Eq. (13) due presence of the black hole are in detail in the 
subsection below. 

The maximum flux will be clearly in the direction of the galactic center. 
The relevant quantity for a measurement is, rather than J(0), the integral 
of J (ip) over the solid angle given by the angular acceptance of a detector 
which is pointing towards the galactic center. We consider therefore the 
function 

(j(o))(AJ7) = -i- / j(v^) dn (14) 

Ai2 J An 

where AO is the angular acceptance of the detector. Fixing for an ACT 
detector a typical value of AQ = 10^^, in Fig 
the core radius a for the allowed regions in Fig 
the galactocentric distance. As can be seen the NFW profile can give very 
high values of ( J (0) ) compared to the other halo models. In Fig. ^ we 
plot instead ( J (0) ) as a function of AQ for the conservative halo models 
shown in Fig. ^ together with the models corresponding to the maximum 
and minimum core radius allowed in Fig. |5| for each profile. 

3.1 Effects of a central black hole 

A small complication occurs when we try to compute the line-of-sight inte- 
gral directly towards the galactic center. There are strong indications that 
a black hole of mass ~ 10*^ Mq resides near the center of the Milky Way. 
The interplay between this black hole and the dark matter halo is difficult 
to model - certainly it is beyond the scope of this paper. One effect of the 
black hole, however, is to gravitationally lens 7 rays that originate from be- 
hind the galactic center, in a cone of angular size on the order of that of the 
Einstein ring. This is given by (see, e.g., [p9| ]) 



UGM Dds .... 
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where G is Newton's constant, M the mass of the black hole, Dds-, D^, Dg 
are the lens-to-source, lens and source distances, respectively. 

In general, for a point source, a double image is produced which usu- 
ally cannot be resolved (the situation is similar to that of microlensing of 
MACHOs in the halo). The resulting magnification may, however, be de- 
tectable and is given by 



+ 2 

+ 4' 



(16) 



where u = P/Oe-, with (3 the undeflected angle to the source. Since gravi- 
tational lensing is a purely geometrical effect, it is of course also operative 
for 7 rays. Returning to our line-of-sight integral Eq. (|l3|), we can now for 
angles if) less than the Einstein angle 9e{1) = do\/{l — Rq)/1 (with around 
1 arcsec for M = 10^ and i?o = 8.5 kpc) split the integral into two parts, 
from to the distance of the galactic center Rq and from Rq to the end of 
the halo. The second part contains the amplification factor /i(Z), thus for 
small angles we replace J{ip) by K{ijj) where 



1 



1 



.5 kpc V0.3GeV/cm3 



Ro 



'Ro 



Oe{1) 



+ 2 



p\i) di{ij) + 



Oe 



(17) 



In Fig. ^ this expression has been used for the smallest angles displayed. As 
can be seen, it causes a mild enhancement for the most non-singular halo 
models and cannot compensate the turnoff of the rise of the NFW profile 
due to the black hole itself. 



4 Present and Planned ACT Detectors 

The present high energy gamma-ray experiments, the Energetic Gamma- 
Ray Experiment Telescope (EGRET) and the Whipple 10m atmospheric 
Cherenkov Telescope (e.g. [^]), lack the sensitivity to detect the annihilation 
lines fluxes predicted for most of the allowed super symmetric models and 
halo profiles. Our aim is to show that the order of magnitude improvement 
in the fiux sensitivity and the substantial improvements in angular resolu- 
tion in the next generation of both ground based and satellite gamma-ray 
experiments, will allow to explore large portions of the MSSM parameter 
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Table 2: Characteristics of Atmosferic Cherenkov Telescopes^ 





Whipple 


GRANITE-III 


VERITAS or HESS 


Effective Area (cm^) 


S-SxlO"* 


5x10'' 


7x10" 


Energy Resolution [oe/E) 


30% 


20%" 


15%' 


Angular Resolution [ag] 


0.14 


0.1^= 


O.O?'' 


Energy Threshold 


250 GeV 


150 GeV 


50 GeV 


Field of View (sr) 


0.001 


0.004 


0.004 



" Estimated from Ref. 
' Resolution at 100 Ge\r 
" Estimated from CAT observations of the CRAB nebula ||34|] . 

Angular resolution at 300 GeV from simulations of a 4 telescope array 
^ Average quantities estimated for the energy ranges: E-y >300 GeV for Whipple; 
E^ > 150 GeV for GRANITE-III; E^ > 75 GeV for VERITAS or HESS 



space searching for the annihilation-line flux, if a halo model which pro- 
vides an increase of the dark matter density towards the galactic center is 
considered. 

In this section we focus on Atmospheric Cherenkov Telescopes (ACTs), 
such as an upgraded version of the Whipple detector currently under con- 
struction (GRANITE-III) as well as proposed arrays of ground-based atmo- 
spheric Cherenkov detectors VERITAS [^] and HESS |33]. We will consider 
in Section |^ the potentiality of the proposed Gamma-Ray Large Area Space 
Telescope (GLAST) experiment. Ground based and satellite experiments 
are to some extent complementary for our purpose. While satellite exper- 
iments offer a lower energy threshold, superior energy resolution, and the 
possibility of relatively long exposure times, the relatively small effective 
area ~ 1 m^ limits the sensitivity at high energies. On the other hand 
ACTs offer the potential for very large effective areas ~0.3 km^, but at a 
higher energy threshold {E^^ > 250 GeV). 

We summarise in Table § the approximate characteristics for existing 
and proposed ACT gamma-ray observatories. 

For the observation of the galactic centre which is treated as a point 
source, we compute the minimum detectable flux N.y{E) using the prescrip- 
tion that, for an exposure of t seconds made with an instrument of effective 
area A^g and angular resolution ag (corresponding to a 68% acceptance of 
events) with corresponding angular acceptance Ail = vrcig, the conditions 
for detection are that the significance of the detection exceed 5ct and that 
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the number of detected photons exceeds 25 events: 

Significance = ^^^^"^^ ^ > 5 (18) 

N^Aest > 25. (19) 

Here, diVbg/dO is the integral number of background events (per unit solid 
angle) falling under the annihilation line. If the dominant source of back- 
ground has a differential spectrum cPN^^g/dEd^l = NqE"^ , and if the energy 
resolution of the instrument is aE/E, then the background under a line at 
energy Eq (i.e., in the interval [Eq — aE,EQ + as] containing 68% of the 
signal) is given by 

where 



ri{aE/E, 6) = \{1 - aE/E)-'+^ - (1 + aE/E)-'+^ 



(211 



gives the the background reduction relative to the integral background spec- 
trum (e.g., r] = 0.35 for an energy resolution of 10%). 

Atmospheric Cherenkov detectors employ large (~10 m) optical reflec- 
tors to image flashes of Cherenkov light from electromagnetic showers formed 
as high energy 7-rays interact in the earth's atmosphere. The Cherenkov 
light pool from a 7-ray (at normal incidence) covers an area of roughly 
5x10^ cm^ on the ground, and defines the sensitive area of the detector. By 
making use of the distinctive differences in the shower images, a 7-ray signal 
can be extracted from the large isotropic background of hadronic showers 
(e.g. [^]). The point of origin of each shower can also be uniquely recon- 
structed from the orientation and shape of the shower image to a precision 
oi a K, 0.14° for a single telescope (e.g. [^]), or roughly a factor of \/iV 
better when N telescopes provide stereoscopic views of the same shower im- 
age. The total level of detected Cherenkov light is roughly proportional to 
the energy of the primary gamma-ray shower and provides an energy reso- 
lution typically of (Te/E ~ 20-40%. A further improvement in this energy 
resolution may be realized with stereoscopic measurements which determine 
the height of the maximum shower development. Typically observations are 
made in a differential mode where each observation on-source is followed by 
an observation of equal duration displaced in right ascension from the source 
direction, so that the control observations sample the same range of azimuth 
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and elevation angles. The number of candidate gamma-ray events is taken 
to be the difference between the number of events passing the data-selection 
criteria on-source and in the control observations. Using this technique one 
can determine the two-dimensional distribution of gamma-ray events on the 
sky, and can reconstruct the energy spectrum of these events assuming that 
the control region does not contain a significant contribution from diffuse 
gamma-rays or from a number of point sources. 

The energy threshold of atmospheric Cherenkov detectors is determined 
by the condition that the Cherenkov light signal must exceed the level of 
fluctuations in the night sky background light. Since both the signal and the 
background light level are proportional to the mirror area A^, the signal to 
noise ratio depends on ^/A^ and the energy threshold is inversely propor- 
tional to this quantity. Since the duration of the Cherenkov pulses is very 
small (on the order of a few nsec), an additional reduction in the energy 
threshold can be realized by decreasing the signal integration time. For the 
largest imaging telescope (the Whipple 10m reflector) the energy threshold 
is approximately 250 GeV. The use of arrays of telescopes operated in co- 
incidence will result in a larger effective mirror area which when combined 
with the use of faster electronics and higher resolution cameras will result 
in a reduction in the energy threshold to perhaps 50 GeV for an array of 
10m telescope. 

Arrays of telescopes also offer the potential for increasing the effective 
area and improving the background rejection. N widely spaced telescopes 
operated independently result in an increase in the effective area by A'^. As 
the telescopes are moved closer together, the effective area is reduced but 
the hadronic rejection can be dramatically improved through better char- 
acterization of the shower development with stereoscopic imaging. For the 
VERITAS (HESS) arrays, 9 (16) 10 m telescopes will be used in a configu- 
ration where each gamma-ray event within the sensitive area will result in 
good shower images in four of the telescopes in the array. In such a mode of 
operation, we estimate a total effective area of 7 x 10^ cm^ and an additional 
hadronic rejection by a factor of ~ 4 through improvements in angular res- 
olution and by another factor of ~ 4 through better characterization of the 
shape of the gamma-ray shower (e.g. Ref. []33| ). 

The sensitivity of atmospheric Cherenkov detectors is determined by 
a relatively large background of misidentified gamma-like hadronic show- 
ers and from cosmic ray electrons, which dominate at lower energies. The 
contribution of the diffuse gamma-ray background is included in the sen- 
sitivity calculation, but has a relatively steep spectrum E^"^-"^ and does 
not contribute significantly to the background for atmospheric Cherenkov 
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telescopes. Even if the diffuse spectrum of the inner galaxy is significantly 
harder than the cosmic-ray spectrum |3^, this conclusion is unchanged. 

Prom data taken with the Whipple 10 m telescope, the measured event 
rates (after rejection of background) and the effective angular aperture and 
effective area (determined from Monte Carlo simulations) can be combined 
to derive the background rate for gamma-like hadronic showers: 



had 



(E >Eo)=6.1x lO-^e^ad 



En 



1 GeV 



-1.7 



cm sec sr 



(22) 



The factor of ehad is introduced to account for improved hadronic rejec- 
tion in a future detector compared with the current Whipple detector. As 
the energy threshold is reduced, the relative efficiency of hadronic showers 
for producing Cherenkov light is reduced with respect to gamma-ray show- 
ers. While this effect would improve the sensitivity at low energies, the 
background rejection from the image analysis generally deteriorates. For 
simplicity we have neglected these effects. 

The showers initiated by cosmic-ray electrons are indistinguishable from 
gamma-rays, and these events can only be rejected on the basis of their 
arrival directions. Cosmic-ray electrons have a steeper spectrum than that 
of cosmic ray nuclei and become the dominant background at lower energies. 



From Ref. |38] the integral spectrum for electrons is 



dn 



-{E > Eo) = 3.0 X 10" 



1 GeV 



-2.3 



cm ^sec ^sr ^ 



(23) 



(at E = 250 GeV, the electron flux is a factor of 5 below the hadronic 
background.) 

Combining the hadron and electron background and taking into account 
the small contribution due to the diffuse gamma ray background, the total 
background under the annihilation line at Eq, is 



diVbg 

dn 



X An = 8.5 X 10^7? 



E, 







Ei 







100 GeV 



-2.3 



(ehad + 0.033) 
3.5 X 108cm2 / V lO^sec 



-1.7 



100 GeV 

t \ r An \ 



+ 0.32- 



O.OOlsr 



) (24) 



This can be used together with Eqs. 18 and 19 to derive the flux sensitivity 
curves of the Whipple atmospheric Cherenkov detector as well as the other 
proposed detectors, as shown in Figs. U and IC. 



The Whipple detector is located in the northern hemisphere, and the 
most likely sites for the VERITAS and HESS arrays are also in the northern 
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hemisphere (although a southern hemisphere site for the HESS array is being 
considered). Since the galactic center is a southern source, it transits at large 
zenith angles for a northern hemisphere observatory (e.g., around 61° at the 
Whipple site). For such observations at large zenith angles, the point of 
maximum shower development is more distant than it would be for showers 
with normal incidence, and the Cherenkov light pool spreads out over a 
larger area. This has the advantage of increasing the effective area but the 
disadvantage of raising the energy threshold. For observations of the galactic 
center at the Whipple site, this results in an increase in the energy threshold 
and effective area by roughly a factor of 5 [39|. For this reason we do not 
consider here the CAT telescope which is located at a latitude of 42° N [^] 
or the HEGRA array ||4^ which will eventually achieve an energy threshold 
of 500 GeV at the zenith but a threshold of 2.5 TeV for the galactic center, 
since in both cases their energy thresholds are too high to be sensitive to 
the majority of parameter space. The 3.8m CANGAROO telescope in the 
southern hemisphere [41| has an energy threshold of approximately 1 TeV 
and will also not probe much of the predicted parameter space. However, the 
CANGAROO collaboration is currently constructing a 10 m telescope which 
may eventually be part of an array of telescopes. This would prove to be one 
of the most important instruments for galactic center observations. Since the 
details of this proposal are not well established, we consider here a generic 
"southern array" which has roughly the characteristics of VERITAS. Since a 
source spends the majority of its time near transit, for such a southern array 
we do not include the effect of large zenith angle observations on increasing 
the effective area. 



5 Prospectives for the Detection of the Monocro- 
matic 7 Line with an ACT 

As we have already pointed out, the detection prospectives for the monocro- 
matic 7 ray line from neutralino annihilation in the galactic halo critically 
depend on the profile which describes the dark matter distribution. This 
is in particular evident for detectors like ACTs which have a small angular 
acceptance. If a large fraction of the total flux emitted is concentrated in 
a tiny region of the sky, the galactic center, whose coordinates are known 
with sufficient accuracy, an ACT can be the ideal instrument for detecting 
neutralino dark matter. This is true for singular halo profiles. We will show 
that with the NFW profile, which is not the most singular profile among 
those that are currently being studied, a signal might be seen already in the 
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next generation of ACT detectors. 

In Fig. P (a) and (b) we show the results of the 27 and Z7 hne fluxes 
for the supersymmetric models considered in Section |2| and for the NFW 
profile which gives the maximal value of ( J (0) ) (10~^ sr) in Fig. We have 
chosen such a small solid angle, which corresponds to the minimal angular 
acceptance for the observation of a point source for an ACT, because for 
a singular halo profile the ratio signal to square root of the background 
increases going to smaller angular acceptances. In the figures we draw also 
the sensitivity curves for each of the ACT detectors considered in preavious 
section. As can be seen, a significant number of supersymmetric models 
could give a signal which exceeds the sensitivity of the Southern Array and 
the VERITAS or HESS detector. The detector in the southern emisphere, 
having a low threshold (around 50 GeV), explores a region of the MSSM 
parameter space which may be accessible to other detection methods, in 
particular to future accelerator experiments. VERITAS or HESS have a 
higher sensitivity for heavier neutralino masses, in an interval which will be 
hardly accessible to direct or other indirect detection methods. 

As the energy resolution of the typical ACT is not better than around 
10%, in the high mass rage the two monochcromatic 27 and Z7 lines, re- 
spectively at the energy ~ and ~ M^{1 - My^M^), are not 
resolvable. In Fig. 10 we sum the contributions to the flux from the two 
annihilation processes. As we have pointed out in Section |2|, a very interest- 
ing feature is that due to Zj annihilation line, the values for the predicted 
flux are concentrated in a quite narrow band and this clearly enhances the 
possibility of a discovery. 



6 The GLAST satellite detector 

In the existing EGRET, and the planned GLAST satellite experiment, high 
energy gamma-rays interact to form an electron positron pair in a distributed 
converter/tracker layer. The measured direction of electron and positron 
in the pair (which improves with increasing energy) and the subsequent 
electromagnetic cascade in the calorimeter layer give the arrival direction 
and energy of each gamma-ray. An anticoincidence shield effectively rejects 
background cosmic-ray nuclei and electrons, and the sensitivity of these 
detectors is limited only by counting statistics and by the diffuse gamma- 
ray background. 

For an energy resolution of 4%, an angular resolution of 0.1° (at E^^ > 
10 GeV) it is clear that the limit given in Eqs. ([l^) and (|T9| ) applies for the 
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point source sensitivity, and that such observations are counting-statistics 
limited. Given this fact it is advantageous to change the observing strategy 
and to relax the angular acceptance window in order that the number of 
collected gamma-rays is increased above the threshold for detection (around 
10 events). 

The estimate of the diffuse gamma ray background is then crucial to give 
a prediction for a signal in a satellite based experiment (for ground-based 
experiments the background from misidentified protons is generally higher). 
The diffuse emission of gamma rays in the galaxy is thought to be due 
mainly to cosmic-ray protons and electrons interacting with the interstellar 
medium, and has an overall enhancement towards the inner galaxy and the 
galactic disk. Given that the galaxy is essentially transparent to gamma 
rays up to energies of about 100 TeV, it is possible to derive the gamma ray 
spectrum in the energy region of our interest from interstellar mass models 
of the galaxy |42, U]. This is however beyond the scope of the paper, we 



will limit ourselves to an order of magnitude estimate which can be obtained 
extrapolating from existing data at lower energies. 

The EGRET experiment has mapped the diffuse gamma ray emission 
up to about 10 GeV. We assume a power law fall-off in energy in the form: 



where we have supposed that the normalization factor A'^o depends only on 
the interstellar matter distribution and is a function of the galactic coor- 
dinates (/,6) only. A^'o has been fixed using the EGRET data at 1 GeV 
published in Ref. [^], in the simple functional form 

Noil,b) =-. r^=S5J 5 + 0.5 if|l|>30° 

= , ^5.5 ^ 0.5 if |1| < 30° (26) 

where the longitude / is assumed to vary in the intervall —180° < / < 180° 
and the latitude b in —90° < b < 90°. The parametrization is in resonable 
agreement with the data at least in the region towards the galactic centre, 
in which we will use it. The slope parameter a may in principle depend on 
both energy and the galactic coordinates. We will, however, adopt a = —2.7 
in lack of data for the energies of interest to us. 

6.1 Acceptance and energy resolution 

The proposed GLAST satellite detector will have a much smaller area (around 
1 m?) than an ACT. However, the solid angle coverage will be of the order 
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of one steradian and the energy resolution will be excellent, at the per cent 
level for gamma rays that enter on the side of the detector and thus pass 
many radiation lengths of the electromagnetic calorimeter. 

As an example [^], the sensitivity of one version of the GLAST calorime- 
ter, consisting of 25 towers (each 32 by 32 cm^) of 10 radition lengths of 
Csl, has been simulated using the GEANT program. The geometric ac- 
ceptance of the instrument was calculated for all gamma rays which fulfil 
that a cylinder around the electromagnetic shower of more than 18 radi- 
ation lengths and a radius of at least two Moliere radii is fully contained 
within the instrument, giving an energy resolution of better than 1.5% for 
gamma rays with energies higher than 50 GeV. It was assumed that the 
diffuse gamma rays enter the instrument from any direction in the upper 
hemisphere, (as the lower hemisphere will be completely screened by the 
Earth). This results in a geometric acceptance of 2.1 m^sr. 

Taking into account the screening by the Earth of the galactic center 
region half of the orbiting time, the useful geometric acceptance in a 1 sr cone 
towards the galactic center is 0.18 m^sr. This is the geometric acceptance 
both for a monoenergetic neutralino annihilation signal from a region around 
the center of the galaxy as well as the diffuse gamma ray background, since 
background photons arriving from other directions can be discriminated by 
the direction of incidence. 

A similar analysis has to our knowledge not been performed for the al- 
ternative design SIFTER, but it appears that there the somewhat larger 
effective area for a given weight of the detector can give a similar improve- 
ment of the 7 ray line discovery potential |4(;]. 



We display in Fig. 11 (a) and (b) the results of the 27 and Z7 line flux 
from Section |2| using the GLAST parameters just discussed, and assuming 
a two-year exposure. Also shown is the curve giving the minimum number 
of events needed to observe an effect at the 5a level, using Eq. ( 25) for the 
galactic diffuse 7-ray background flux. The NEW halo model corresponding 
to the maximal rate was assumed for the halo. 

As can be seen, a fair fraction of our set of MSSM points can be probed 
under these circumstances. An advantage with a detector of large acceptance 
is that "clumps" of dark matter of relatively small angular size could still 
be detected. Since the annihilation flux is proportional to the square of the 
dark matter density, such enhancements would appear as bright spots on the 
7 ray sky, without optical counterpart. Another advantage with a satellite 
experiment is the excellent energy resolution possible, which may enable 
a separation of the two lines from 27 and Z7 respectively. If both lines 
were to be observed, a comparison of line strengths would give interesting 
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information on the super symmetric model. 

If the relative energy resolution is e = aE/E, the two lines are separable 
provided that 

My 

M,S-|, (27) 



which for e = 0.015 gives < 370 GeV. 



7 Continuum Gamma Rays 

Besides the gamma ray line coming from the 77 and Z7 final states, there 
may also be produced a continuum 7 spectrum mainly from decaying tt^ 



mesons created in the fragmentation of quarks 49 1. The general 

drawback of the continuum 7 spectrum is that in contrast to the 7 lines 
it lacks distinctive features which can exclude other more mundane sources 
of production if a signal is believed to be found. On the other hand, the 
number of photons from this continuum source can be much higher than 
from the line processes. 

In fact, should a 7 ray line be found, by necessity there has to be a con- 
tinuum signal also, although in some cases severely affected by the galactic 
diffuse background. We therefore investigate what can be expected for the 
higgsino-like high-mass sample, since previous treatments were mainly done 
for low-mass neutralinos. 

For our high-mass sample, the main annihilation mode is into a WW or 
ZZ pair. We have simulated the continuum photon spectrum from these and 
qq final states through the PYTHIA 6.113 Monte Carlo code ||5^. Introduc- 
ing the scaling variable x = E^/M^, we find empirically that the quantity 
x^'^dN^/dx is well described by an exponential in the variable x, as shown 
in Fig. H for = 500 GeV. 

As can be seen, t and b quarks give a quite soft spectrum. Since u quarks 
(as well as s and d quarks) have a very small branching ratio for these high 
mass neutralinos, the dominant source of photons above x ~ 0.1 will be WW 
and ZZ. Both of these final states show a very simple scaling behavior when 
increasing the energy (i.e. the neutralino mass). This is clearly shown in 



Fig. |13|, where photons from both WW and ZZ at of both 500 and 2000 
GeV can be described to the level of accuracy required for our purpose by 
the distribution 

dN^ _ M^dN^ _ 7.8, 

dx ~ dE^ ~ x^-^ ■ ^ ' 
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The models in our sample with > 300 GeV with the highest rate of 
WW and ZZ production have a rate given approximately by 

- 5 . 10-2^ (^^^)' cm^s-i. (29) 
The observed differential flux at the Earth is then given by 



'300 GeY V dN^ 



■J{^!) cm^^s^^sr^^GeV"^ , (30) 



where J{^) is defined in Eq. (p^). 

Using the NEW halo profile which gives the highest rate, and integrating 
over Ari = 10~^ sr towards the galactic center, this gives the results shown 



in Eig. y for = 300, 500, 1000 and 2000 GeV. The galactic background 



curve is the estimate Eq. ( |25|) , which for 6, / = gives 



~ 6 • 10-^ (^^) cm-2 s-i sr-i GeV-^ . (31) 

We see that below a mass of 1000 GeV, this additional component of 
the photon spectrum is larger than the estimated background over some 
energy range, and could thus be observable, see Eig. ^ where also the 
results for a 50 GeV neutralino (annihilating mainly to hh) pairs are shown, 
using the parametrisation of ||4^ of the 7 spectrum from b quarks. As 
can be seen, at least the 50 GeV candidate would appear to be excluded 
already by EGRET data, were it not for the uncertainties related to the 
halo parameters. The combined features of a break in the spectrum and an 
angular distribution that does not follow that of the disk and bulge of the 
galaxy could make this type of signal distinguishable from the background, 
although it may be difficult to exclude other sources of 7 rays giving similar 
features. In fact, present EGRET observations are not inconsistent with a 
continuum spectrum originating from dark matter annihilations, but other 



explanations are possible as well [51|. 



It is interesting to note in Eig. |15| that a 300 GeV higgsino would cause 
a flatter spectrum than the "canonical" E""^'"^ in the 1 - 100 GeV range. It 
is intriguing that a flatter spectrum is indeed observed towards the galactic 
centre, something which has turned out difficult to reproduce in conven- 
tional models [^]. Before detection of a line signal it will, however, not be 
possible to rule out other sources of high energy 7 rays. In any case, the 
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hitherto unexplored energy band between 100 and 300 GeV will be of crucial 
importance to verify or rule out a supersymmetric component of the 7 ray 
spectrum. 

The higher-mass neutralinos (above 1 TeV) do not give as useful a contin- 
uum signal because of the rapid decrease with of the signal (see Eq. (|30|)), 
whereas the 7 line remains a viable signature up to that region. 

8 Discussion 

To conclude, we have shown that the new high energy cosmic 7-ray de- 
tectors, both ground- and space-based will have an interesting opportunity 
to search for signals of dark matter paticle annihilations in the Milky Way 
halo. Indeed, existing data from EGRET may already rule out some of the 
low-mass MSSM models if the halo is of the most singular form studied here. 

Even neutralinos as heavy as several hundred GeV could give continuum 
7s with an observable flatter spectrum than the E~'^''^ expected in the sim- 
plest models of diffuse galactic 7 rays, in the direction towards the galactic 
center. The absence of a clear signature may, however, make an eventual 
claim of detection of a supersymmetric signal difficult to defend. 

The only reliable 7-ray signals indicating directly the existence of slowly 
moving, heavy neutral particles in the Milky Way halo seem to be the 7-ray 
lines produced by XX ~^ 77 aiid XX ~^ -^7 annihilations. These lines, of 
intrinsic relative width of the order of 10~^, have no plausible astrophysical 
background whatsoever and would constitute a "smoking gun" of supersym- 
metric dark matter. 

As we have shown, the larger rates for the line processes found in the 
recent first full one-loop calculations in the MSSM, together with the central 
enhancement indicated in simulations of cold dark matter halos, make these 
processes accessible to the next generation of 7-ray detectors over a fair range 
of parameters. In particular, if the neutralino is the lightest supersymmetric 
particle and is heavier than around 1 TeV, detection of the 7 line may well 
be the only way to discover supersymmetry in the foreseeable future, as all 
other methods (accelerator, direct detection and indirect neutrino detection) 
are discouragingly far from the required sensitivity. 
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Figure 1: Annihilation rate of neutralinos into the 27 final state for the 
sample of super symmetric models described in Section |2|. A different marker 
color is used for the three classes of neutralinos we have defined: a red marker 
indicates a higgsino-like neutralino, a green marker a mixed neutralino, while 
blue markers are for gaugino-like neutralinos. 
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Figure 2: Annihilation rate of neutralinos into the Z7 final state for the 
sample of supersymmetric models described in Section ^. The marker colors 
have the same meaning as in Fig. ffl. 
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Figure 3: vaz'y versus 2 va2'y for the sample of supersymmetric models de- 
scribed in Section ^. The marker color has the same meaning as in Fig. |l[ 
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Figure 4: 2 1^(727 and vaz-y for a sample of mixed neutralinos in the mass 
range 150 - 400 GeV. For each model the class of diagrams which gives the 
main contribution to the cross section is indicated. 
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Figure 5: Allowed values of the parameters po, the local halo density, and 
a, the core radius, for the four halo profiles considered in Section ^. The 
allowed regions vary as a function of the third parameter which enters the 
discussion, the galacto centric distance of the solar system Rq; we plot the 
regions corresponding to Rq = 8.5 kpc which extend to lower values of a and 
to Rq = 7.1 kpc which allow higher values of a. The markers indicate the 
halo profiles which are considered in Fig. ^. 
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Figure 6: Function J{ip) as defined in Eq. (13) for the four halo profiles 
considered in Section The choice of the parameters a and po is indicated 
in the figure, while we have fixed Rq = 8.5 kpc. 
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Figure 7: ( J (0) ) as defined in Eq. (14) for Ail = lO^^sr and as a function 
of the core radius a for the four halo profiles considered in Section |^ Each 
contour is given by the maximal and minimal values of the parameter pQ as 
given in Fig. |5|, for the two choices of Rq. 
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Figure 8: ( J (0) ) as a function of AO for the halo profiles defined by the 
maximal and minimal values of the core radius allowed in Fig. ^ and for 
those considered in Fig. ^. Curves for both choices of Rq are shown. 
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Figure 9: Gamma ray flux from a 10~^ sr cone encompassing tlie galactic 
center for the 27 (on the left) and the annihilation line (on the right). 
The NFW halo profile giving the maximal flux has been assumed. The solid 
lines show the 5a sensitivity curves of the ACT detectors described in the 
text. 
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Figure 10: Gamma ray flux from a 10~^ sr cone encompassing the galactic 
center summing the contributions of the 27 and the annihilation lines 
for heavy neutralinos. The NFW halo profile giving the maximal flux has 
been assumed. 
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Figure 11: The number of events expected in GLAST from a 1 sr cone en- 
compassing the galactic center, assuming a 2 year exposure and calorimetry 
as described in the text, for the 27 (on the left) and the Z7 annihilation 
line (on the right). The NFW halo profile giving the maximal flux has been 
assumed. The solid line shows the number of events needed to obtain a 5a 
detection over the background as estimated from EGRET data. 



38 



0.1 



0.01 





1 


1 1 




1 








Continuum Photons from 500 GeV 






+ 




Neutralino Annihilations 




































• 


u quarks 




" . + • 




■ 


W bosons 




•. » 




▲ 


Z bosons 




■ i. '4 
■» + 







b quarks 






+ 


t quarks 






■ • 


• 










+ 1 • 

■ 


• • 








- 


* ■ 

+ 


• • 

1 • • • 










+ 

+ 


i 4 ■ • . 
i ■ • . 














• 

• • 


• 








+ ( ( i ■ , 


• 

• 


• 








J- ' * ■ 


■ 












■ 

■ ■ 

A ■ 


• 

• 

• 








6 


* A 


■ ■ • 

■ 


• 






+ * 


A A 










+ 


« 




• 






+ 

+ 






■ - 
■ 






+ 

+ 


+ 


( « 
1 « 


A 

A 


4- 


1 


1 1 




1 


4 



0.0 



0.1 



0.2 0.3 

X 



0.4 



0.5 



Figure 12: Continuum photons from 500 GeV neutralino annihilations. The 
rate is given as x^^'^dN^ / dx per final state pair for u quarks (filled circles), W 
bosons (filled squares) , Z bosons (filled triangles) , h quarks (open diamonds) 
and t quarks (crosses). The scaling variable x is defined as a; = E^/m^, 
where E^/ is the photon energy and 2m^ the total energy of the final state in 
the annihilation. Note that for a realistic neutralino of this high mass, the 
annihilation into u quarks, normalized to unity in the figure, will generally 
be much suppressed. 
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Photons from W and Z bosons in 
Neutralino Annihilations 
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Figure 13: Continuum photons from W and Z bosons in in high-mass neu- 
trahno (in practice, mostly higgsino) annihilations. Pythia [50| results for 
500 GeV neutralinos are shown with filled symbols, for 2000 GeV with open 
symbols. The line gives the analytic fit according to Eq. (28). 
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Figure 14: Continuum photons from the direction of the galactic centre 
originating from W and Z bosons in neutralino annihilations, for = 
300, 500, 1000 and 2000 GeV. The NFW halo profile giving the maximal 
signal has been assumed, and an angular integration over 10"'^ sr performed. 
The background flux (solid line) is that predicted by Eq. (^). 
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Figure 15: Total photon spectrum from the direction of the galactic centre 
originating from W and Z bosons in neutralino annihilations, for 300 GeV, 
and from b quarks for 50 GeV. The NFW halo profile giving the maximal 
signal has been assumed, and an angular integration over 10~^ sr performed. 
The background flux is that predicted by Eq. (|25|). 

In addition, the maximal 7 line strength found in our sample is displayed 
for these two masses, assuming a relative line width of 10~^. 
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